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Abstract. In this paper we show normalized differential source counts n(S) at 408 MHz and 1.4 GHz of radio 
sources separately for FRI and FRIT classes with extended and compact morphologies. The maps from the FIRST, 
NVSS, and WENSS surveys are used to define the source morphology and flux density. The counts provide a basis 
for a direct test as well as constraining the cosmological evolution of powerful extragalactic radio sources in terms 
of the dual-population model (Jackson & Wall 1999), where radio sources of Fanaroff-Riley (1974) types I and 
II are regarded as two physically separate types of active galactic nuclei (AGN). The predicted count values are 
compared with the observational data to find the best fits for the evolution and beaming parameters, and to 
further refine the model. 
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1. Introduction 

Galaxies (and/or quasars) hosting active galactic nuclei 
(AGN) are powerful extragalactic radio sources which pro- 
duce jets and extended radio emitting regions (lobes) of 
plasma (Scheuer 1974). Most of these sources can be easily 
separated into two distinct classes: (i) the edge-brightened 
FRII sources with well-collimated jets, radio lobes and 
prominent hot spots, and (ii) edge-dimmed FRI sources 
with bending jets and peaked radio emission near the nu- 
cleus. There is a clear division in the luminosity between 
FRI and FRII, which lie, respectively, below and above 
Li78 « 2x10^5 WHz-^sr-i (Fanaroff & Riley 1974). Owen 
& Ledlow (1994) demonstrated that this separation is re- 
produced in the two-dimensional radio-optical luminosity 
plane in an unambiguous way. Furthermore, it is now quite 
clear that FRII and FRI sources differ radically their cos- 
mological evolution (Wall 1980), i.e. the density and/or 
luminosity evolution rate of FRII is high, while that of 
FRI seems to be rather low. 

Radio source counts n(S) involving distributions of the 
other observational parameters (i.e. luminosity, spectral 
index, redshift) were used to constrain a number of nu- 
merical models of cosmological evolution of extragalactic 
radio sources (e.g. Wall et al. 1980; Peacock & Gull 1981; 
Condon 1984; Peacock 1985; Dunlop & Peacock 1990). All 
these models assumed two distinct populations of sources: 
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steep-spectrum and flat-spectrum ones, which in principle 
were selected based on their spectral characteristics only. 
Since the late 70s, however, there have been efforts to 
find a common scheme for some types of extragalactic ra- 
dio sources (Rowan-Robinson 1977; Scheuer & Readhead 
1979). These unification attempts aimed to explain the 
diversity of observational properties of extragalactic ra- 
dio sources in terms of one simple model. The starting 
point for such models is the intrinsic structure of a radio 
source. Because they are not spherically symmetric, their 
observed form depends on the angle of view. The realiza- 
tion that orientation angle in AGN is of great importance 
for unification came with the discovery of observational su- 
perluminal motion in some compact radio sources (Cohen 
et al. 1971; Moffet et al. 1972) and its interpretation in 
terms of bulk relativistic movements of the jet's plasma 
(de Young 1972). It is obvious that when some emitting 
body moves in a relativistic manner, the radiation received 
by an observer is a function of the angle between the line 
of sight and the direction of motion. The small fraction of 
objects with compact morphologies are those viewed with 
their axes inclined at small angles to the line of sight. 
Furthermore, the beamed sources should be linked to an 
unbeamed parent population. The unified scheme of AGN 
has been quite successful in explaining differences between 
various classes of radio galaxies on the basis of the jet lu- 
minosity, its speed, relativistic beaming effect, orientation 
towards the observer, and environmental effects (for a re- 
view, see Urry & Padovani 1995). 
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Recently, Wall and Jackson (Wall & Jackson 1997; 
Jackson & Wall 1999) [hereafter WJ] introduced the dual- 
population unified scheme. Their model consists of two 
parent populations of powerful extragalactic radio sources, 
FRI and FRII, which appear (dependent on the viewing 
angle) as lobe-dominated (LD) steep-spectrum sources, or 
core-dominated (CD) flat-spectrum sources. 

Our main goal is to test the dual-population model 
of FRI and FRII radio source cosmological evolution and 
to refine its parameters with new well-defined samples of 
sources at two radio frequencies and within a wide range of 
radio flux density. The recent radio surveys with sufficient 
sensitivity, spatial resolution and fltix-density limit, i.e. 
FIRST (Becker et al. 1995), NVSS (Condon et al. 1998) 
and WENSS (RengeHnk et al. 1997), allowed us to deter- 
mine morphologies, fltix densities, and spectral indices for 
a large enough number of objects and thereby proved to 
be suitable for constructing such samples. 

In the following section the complete source samples 
used in this analysis and the counts of extragalactic ra- 
dio sources n(S) based on data at 1.4 GHz and 408 MHz 
are described (Sect. 2). In Sect. 3 the application of the 
dual-population numerical model of radio source evolu- 
tion is shown. In Sect. 4 the results of modelling of the 
counts n(S) compared to the observational data and dis- 
cussion arc given. Finally, Sect. 5 contains the summary. 
Throughout the paper we assume Hq = 50 kms~^Mpc~"'^ 
and the Einstein-de Sitter cosmological model (Q = fim = 

!)• 

2. Samples and observational source counts 

To construct observational source counts we selected 
two complete samples of extragalactic radio sources. The 
1.4 GHz sample is selected from four different radio sur- 
veys, namely the FIRST and Mitchell (1983) surveys, the 
GB/GB2 (Machalski 1998), and BDFL (Bridle et al. 1972; 
Bridle fc Fomalont 1974) samples. The 408 MHz sample 
is based on the third Bologna radio survey (Ficarra ct al. 
1985). We exclude all starburst as well as galactic sources 
from our study. The detailed description of the 1.4 GHz 
and 408 MHz samples as well as the criterion of the source 
selection are described in Sect. 2.1 and Sect. 2.2, respec- 
tively. 

We divide all sources from our two samples into four 
distinct types, i.e. (i) lobe-dominated edge-dimmed ob- 
jects (LD-FRI), (ii) lobe-dominated edge-brightened ob- 
jects (LD-FRII), (iii) compact objects with steep spectra 
(CSS), and (iv) compact objects with flat spectra (CD). 
A compact object is counted as CD-type when its low- 
frequency spectral index a (if S^ oc where S^ is 
the flux density at frequency v) is less than 0.6 (details 
about spectral index determinations are given below). Our 
CSS-class of radio sources includes also GigaHertz-Peakcd 
Spectrum (GPS) sources which appear in the samples. In 
addition, in our further analysis we treat the CSS-type ob- 
jects together with the lobe-dominated FRII ones, because 
it is believed that they are just young counterparts of the 



older classical doubles (de Vires et al. 2002). Furthermore, 
as there is no room for lobe-dominated transient morphol- 
ogy (FRII/FRI) objects in the WJ scheme, we put all the 
"hybrid" sources into the LD-FRII class. In general, in 
this simple approach we have only three groups of radio 
sources: LD-FRI, LD-FRII and CD sources. The number 
of sources of a particular type within the selected 1.4 GHz 
and 408 MHz flux-density ranges is given in Table 1 and 
Table 2, respectively. 

2.1. 1.4 GHz sample 

The 1.4 GHz sample consists of objects with observed flux 
densities between 20 mJy and 20 Jy (the flux scale of Baars 
et al. 1977 is adopted) which, depending on their flux, were 
selected from different surveys. 

Weak source subsample. 

The low flux-density subsample consists of sources with 
20 < Si. 4 < 200 mJy, selected from the FIRST and 
Mitchell (1983) radio surveys. Both surveys were made 
using the Very Large Array with 5" and 15" angular reso- 
lution respectively. The flux density of each object is taken 
from the NVSS source list. To determine morphological 
types of the extended radio sources we create a contour 
map for each of them using the FIRST and/or NVSS sur- 
veys. Next we assign by eye the particular objects in the 
maps to FRI or FRII type. The completeness of the FIRST 
and NVSS surveys is 100% in the above flux-density range, 
however, because of their limited baselines they are insen- 
sitive to extended structures (i.e. larger than 15' in the 
case of NVSS). To differentiate between compact steep- 
spectrum (CSS) and flat-spectrum (CD) sources we cal- 
culate a spectral index between 325 MHz (WENSS) and 
1.4 GHz (NVSS). In some cases when the source is not 
included in the WENSS catalog, its low-frequency flux 
density is taken from other radio catalogs. Because; the 
spectral information about sources comes from the flux 
density measurements performed at two frequencies only, 
some of the compact GigaHertz-Peakcd Spectrum sources 
could be erroneously classified as CD ones. However, the 
fraction of GPS sources at 1.4 GHz constitutes at most 
only about 5% of the entire sample (see Table 1 of Jackson 
& Wafl 1999). 

The sources from the FIRST survey are selected from 

a small region of the sky between 07*^19™ < 02000 < 
OSi^OO" and -h28°45'< ^2000 < +31°10'. The sources from 
the Mitchell survey are selected from a sky region be- 
tween 12'^56'" < ^2000 < 13^20"" and +28°16' < (^2000 < 
+32°02 . There arc 400 radio sources in this subsample. 

Intermediate source subsample. 

This subsample is taken from the GB/GB2 revised list 
whose definition is given by Machalski (1998) and consists 
of 346 sources with fluxes between 200 mJy and 2 Jy. 
All sources have bn > 20°. The value of each source's 
spectral index and the shape of the spectrum is taken 
from Machalski (1998). 
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Table 1. Observational normalized differential source counts at 1.4 GHz. 



AS 

[Jy] 



sky area 



LD-FRI 
log(n(S)) 



LD-FRII 
log(n(S)) 

[sr-ijyi-S] 



CSS 

log(n(S)) 

-ijyi-S] 



sr 



CD 

log(n(S)) 

[sr-ijyi-5] 



0.02 
0.05 
0.10 
0.20 
0.25 
0.55 
1.00 
2.00 
5.00-f 
10.00 



f0.05 
^0.10 
^0.20 
fO.25 
^0.55 
M.OO 
^2.00 
^5.00 
10.00 
^20.00 



0.0133 
0.0133 
0.0133 
0.0550 
0.0906 
0.4414 
0.4414 
4.30 
4.30 
4.30 



25 1.06; 



3 
3 

10 
14 
3 

20 
11 
1 



+0.08 
0.10 
-, 91+0.13 
^•^J^-0.19 
1 99+0.20 
-■-•^^-0.37 
1 4-1+0.20 
^•^^-0.37 

1 S5+°-ii 

^^■■^^-0.13 
1 91+0.20 
-^■^-^-0.38 

i.44t»:?? 

1.24l°-3° 



94 
47 
23 
16 
52 
93 
30 
113 
12 
3 



1 64+°°* 

1 08+0-06 
-'-•^°-0.07 

9 -in+0-09 

2 14+010 
^■^^-0.12 
9 90+0-05 

^•^'J-0.07 
9 oc+0-04 
'^■•^'^-0.05 
9 9-1+0-07 
^•^-■--0.09 
9 jQ+0-04 

i-87i°;lJ 

+0-20 

-0.37 



1.72 



71 

28 
19 
4 
19 
16 
11 
30 
6 
2 



1+0.06 

-0-05 

hO.08 
-0.09 
5+0.09 
-0.11 

1 54+°-!^ 



1.511 

1.751 
2.021 



-0-30 

-,+0.09 
-0.11 
hO.lO 
-0.12 
7+0.11 
-0.15 

1 62+°°'' 

-■-•"^-0.09 
1 c;7+0.15 
-■-•^'-0.23 
1 ci4+0.24 



1.79: 

1.581 
1.771 



60 

13 



27 
27 
13 
42 
5 
2 



1 44+0-05 

^■^^-0.06 

1 42+°" 

1-^^-0.14 

1- "-'-0.18 
1 S/l+0-13 
^■^'^-0.19 
1 94+0-08 
-'-•^^-0.09 
, +0.08 



0.09 



I.8I: 

1 s/i+o-ii 
1 76+°°'' 

1 4Q+0-16 
-■-■^^-0.26 
1 t:4+0-24 



Table 2. Observational normalized differential source counts at 408 MHz. 



AS 

[Jy] 
0.10+-0.20 

0.20^0.40 
0.40^0.80 
0.80^1.60 
1.60^5.00 
5.00^25.00 



sky area n 

[sr] 

0.0102 
0.0327 
0.0785 
0.1773 
0.2774 
0.2774 



LD-FRI 
log(n(S)) 

[sr-ijyi-«] 



12 
10 
9 
9 

5 
2 



1.941 



f 0.11 

-0.15 
-I 01+0.12 
J^-OJ^-0.16 

1 84+°-i2 
1 94+°-^^ 

i.80_Q 25 
9 10+0.23 



LD-FRII 
log(n(S)) 

[sr-ijyi-«] 



62 
91 
95 
81 
58 
10 



9 f.r+0.05 

2.76^°:°^ 
2.86l°;°t 
2 89+° °^ 

^•"^-0.06 

9 00+0-05 

^•■'^-o.os 

+0.12 
-0.17 



2. 



29 
44 
50 
31 
15 
3 



CSS 

Iog(n(S)) 

[sr-ijyi-S] 



2.32I 
2.451 



f 0.07 
^-0.09 
f 0.06 
-0.06 
9 ro+0.06 

9 47+0.07 
-0.09 

9 qo + 0.10 
Z.OO_Q 13 

9 qr + 0.20 
^.ou_Q 3Y 



13 
14 
16 
17 
10 




CD 

Iog(n(S)) 

[sr-ijyi-«] 



07+0-11 

-0.14 

95+0-10 

•'^'-'-0.14 
OQ+0-09 

•"^-0-13 
91+0.10 

•^J^-0-12 
1 c;+0-12 

•l'J-0.18 

< 1.88 



Strong source subsamplc. 

The subsample of 247 sources with flux densities 2 < 
Si. 4 < 20 Jy, found within the area of the sky defined 
by -5° < (^2000 < +70°, bn > 20°, is taken from BDFL 
sample. The original source's flux has been transformed 
from the KPW (Kellermann et al. 1969) to the Baars scale 
(Skpw X 1.029 = Ssaars)- In this subsample the spectral 
index is taken from the Kiihr catalogue (Kiihr et al. 1981). 

2.2. 408 MHz sample 

The 408 MHz sample consists of 686 extragalactic sources 
with flux 0.10 < S408 < 25 Jy selected from the B3-VLA 
catalogue (Vigotti et al. 1989). The original 408 MHz flux 
densities are taken from The Third Bologna Catalogue 
(Ficarra et al. 1985). The sample includes objects within 
the R.A. range 0107*^03" < aigso < IS'^OO™. The spectral 
index for sources from this sample is determined between 
408 MHz (B3) and 1.4 GHz (NVSS). 

2.3. Source counts 

The observational normalized differential source counts 
n(S) are defined as: 



n(S) 



nA- 



ASS-2-5 ' 



(1) 



where n is the number of sources (e.g. of lobe dominated 
FRI or FRII, or of core dominated ones) in the sky area 
A with flux density between Si and S2 (AS = S2 — Si), 
and 



S = exp{(log(Si) + 



log(S2) - log(Si) 



) X In(lO)}. 



(2) 



The errors An(S) are derived as: 
(n±VS)A-i 



An(S) 



ASS-2-5 



n(S)|. 



We also determine the observational fraction F of sources 
of a particular morphology within a given flux-density 
range (AS), 

nAS 

and (0 < F < 1), where uas is the number of all sources 
in the same flux-density range. The error of the fraction 
F is derived as: 



AF 



F(l-F) 



(3) 



The sources from the samples described in Sect. 2 are 
used to construct the separate radio source distributions 
of core-dominated and lobe-dominated FRI or FRII ob- 
jects. The resulting normalized differential counts n(S) 
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at 1.4 GHz and 408 MHz are presented in Table 1 and 
Table 2, respectively. The columns give the adopted flux- 
density intervals, the area of sky, the observed number of 
sources of particular morphology and the relevant differ- 
ential source counts with estimated errors. We are aware 
that our procedure of morphological classification of radio 
sources is not quite perfect and may occasionally lead to 
erroneous assignments. The misclassified cases can have a 
significant effect on the resulting counts when a number of 
objects within the given interval of flux density is small. 
Figs. 1 and 2 show the counts described above, while in 
Fig. 3, the observed source fraction F, both at 1.4 GHz 
and 408 MHz, are plotted. 

3. Dual-population model 

The dual-population paradigm for powerful extragalac- 
tic radio sources, proposed by WJ, involves two parent 
populations: FRII and FRl type sources. The orienta- 
tion angle of the radio axis to the plane of the sky de- 
termines the observed characteristics and classification of 
the objects as lobe-dominated steep-spectrum sources or 
core-dominated flat-spectrum ones. The beaming model 
is applied to those distinct populations to model their re- 
spective contributions to the source counts. The evolution 
and beaming scheme can be used to predict the intensity- 
dependent n(S) distributions for beamed and imbeamed 
counterparts of FRI and FRll radio source populations at 
any frequency. The consecutive steps to model the n(S) 
counts are: 1) constructing the local luminosity function 
and assuming a form for the evolution function, 2) us- 
ing the beaming model to compute the flux-densities and 
spectral indices of modeled sources, 3) computing the ex- 
pected source counts and comparing them to the observa- 
tional data, and 4) varying the values of model parameters 
and repeating the steps above until the predicted counts 
fit well the observations. The above procedure is applied 
to FRI and FRH sources separately. 

3.1. Radio luminosity function and evolution function 

The radio luminosity function p(L, z, z/)dL specifics the 
comoving number density of sources with luminosities L 
to L-|-dL measured at frequency v (in the source frame) 
and which are at redshift z. The local luminosity function 
p(L, zq, ly) along with an evolution function E(L, z), which 
modifies the former one, give the radio luminosity func- 
tion at any epoch p{h, z, = p{h, zq, i^)E(L, z). We adopt 
the form of the local radio luminosity function directly 
from WJ. The evolution function E(L, z) is chosen as an 
exponential luminosity-dependent density evolution, first 
proposed by Wall et al. (1980). The evolution function has 
the form: 

E(L,z) = expM(L)T(z), (4) 
where t is the so-called "look-back" time given by 
t(z) = 1-(1 + z)-i-5. 



We also apply a redshift cutoff Zc, so that for z > Zc, 
E(L,z) = 0. The evolution rate M depends on a source 
radio power L: 



M(L) = M^.. /;«,V-iogL\ 
M(L) = 
M(L) = M„,ax 



3.2. Relativistic beaming 



for Li < L < L2, 
for L < Li, 
for L > L2. 



In the beaming scenario, it is assumed that all radio 

sources are intrinsically stecp-spcctrum, and fiat-spectrum 
sources are merely those in which the core radiation is 
boosted in a sufficiently strong way to dominate the to- 
tal flux density. The beaming model predicts that the 
observed radio flux density S^, from a given source at a 
given frequency v is composed of flat-spectrum emission 
from the core and steep-spectrum from the extended lobes 
emission, S^cd and Si/ld respectively. The observed core- 
to-extended flux ratio R, which depends on the the trans- 
verse (unbeamcd) core dominance parameter, line-of-sight 
orientation of the jet towards the observer and the jet's 
speed, is given by: 



R{9, v) 



Rc(i^){[7(l - /3cos(?)]-P + [7(1 + /3cos0)]-P}, 



(5) 



where Rc(i') is the frequency-dependent intrinsic core-to- 
extended flux ratio, 9 {0° <9 < 90°) is the angle between 
the jet axis and line-of-sight. 7 is the Lorentz factor of the 



plasma, 7 = 



■■. For the radio emission with contin- 



uously ejected plasma, p = 2 — acD, where cccd is the 
spectral index of core emission, for which we assume the 
value of 0.0. Following WJ, we consider that the mean 
spectral index is old = 0.75, and the total flux density of 
a source is S,y = S,yCD + S,yLD = S,yLD(R(6'. ^) + ^)- Tlw- flux 
densities of the flat-spectrum component at two frequen- 
cies vi and V2 are equal (Sj^jCd = S^jCd)- The parameters 
of 7 and Rc(i^) arc the same for each object from a given 
population of radio sources. A source is classified either 
as a core-dominated one (CD) if 9 is small enough so that 
spectral index of the entire source is q;[^J < 0.6, or a lobe- 
dominated one (LD) for aj^^ > 0.6. This spectral index is 
given by the equation (Morisawa & Takahara 1987): 

R(9,i^2)-l-l i 

(6) 



log- 



The critical observed core-to-extended flux ratio 

Rmin(^, i^), dividing the sources into the CD or LD 
type at a frequency 1^2 is given by: 



Rn 



tin\aJl-ai,-D 



1 



1 



\ V2 ' 



(7) 



Rmin(^, v) corresponds to the critical angle 9, which can 
be determined from equation: 

1 



coa{9) 



(3 /37'R> 



l_^RcM^p- 



(8) 
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assuming that the beaming due to the counter-jet is neg- 
ligible. 

3.3. Expected n(S) distributions 

The number of sources with flux-densities S to S-|-dS at 
frequency v are given by: 



n(S, ly) 



z=0 JR 



n(S,z,z/)P(R(6',i'))dRdz, 



(9) 



where n(S,z, j/)dS is the number of sources in An sr of 
the sky with flux-densities S to S+dS at frequency and 
redshift z. P(R) is the distribution of R(0, i/) and in a 
sample of randomly orientated radio sources with the same 
values of Rc(f^) and 7, P{R{9, v)) is given by: 

P(R(6l,iy))dR = d(cos6i), 

so equation ini becomes: 

n(S,j/)=/ / n(S,z,t/)sin6'd6ldz. (10) 

The values of ^min and 0max depend on the source mor- 
phology. For CD sources ^min = 0°, and ^max = 90° for 
LD sources. 6'inax for CD is equal to 9nun for LD sources 
and can be derived from equation |S1 

The radio luminosity function p(L, z, v) and the number 
of sources n(S, z, ly) are related in the following way: 

n(S, z, j/)dS = p(L(S), z, j/)dLdV, (11) 

where dV is a volume element and is equal to 47rD^dr. 

2c{flz z)[{nz + l)i - 1]} 



D = 
and 

dr — 



Ho(l + z)f]2 



(12) 



cdz 



Ho(l + z)(l + f7z)2 



(13) 



Inserting L = D2(1+z)^+"S and eauationlT^into equation 
nil yields: 



n(S,z,:/) = 



47rD'*c(l + z)"p(L, z, z^) 
Ho(l + f7z)^L 



(14) 



It is useful to calculate the normalized counts and for that 
purpose we multiply equation [TUbv S^: 



I 47rc(l + z) 2 2«p(L,z,:/)L2 
Ho(l + f7z)5D 



(15) 



4. Results and discussion 

Following the procedure given in the preceding section we 
derive the normalized differential counts and the source 
fractions. The set of parameter values originally proposed 
by WJ is called 'model A' (see Table 3). It is necessary 
to modify these values to fit the observational data. The 
best- modeled fit to the source counts at 1.4 GHz and 
0.408 GHz is called 'model B' and its parameter values 
are given also in Table 3. The values are determined using 
statistical methods, by minimizing evaluated between 
the observed and modeled source distributions, both at 
1.4 and 0.408 GHz. The value of is evaluated as: 



i=l 



errdata 



where njatai is the value of the data in the ith bin, 
similarly nmodoU and errjatai are the model value and 
data error in the ith bin, respectively. These all are 
determined in the logarithm space and the sum is taken 
over the N bins. 

The observed and calculated n(S) distributions for 1.4 
and 0.408 GHz are presented in Fig. 1 and Fig. 2, re- 
spectively. Fig. 3 shows the expected and observed source 
fraction. 

The best-fitted parameters of cosmological evolution 
and relativistic beaming that we estimate are slightly dif- 
ferent from those obtained by WJ on the basis of the 
5 GHz data. However, it is important to note that our 
calculations do not strictly follow all the steps proposed 
by Jackson & Wall (1999). For example, we do not de- 
mand the evolution function E(L, z) to peak at Zc/2. 
Furthermore, our approach does not involve any depen- 
dence of the intrinsic core-to-extended flux ratio Rc on 
radio power. We estimate that the transition between 
evolving and non-evolving FRII sources appears for lu- 
minosities Li5i = lO^'*'^ (WHz~^sr~^), whereas objects 
with Li5i > 10^'^-^^ (WHz~^sr~^) show maximal evolu- 



tion with the rate of Mj, 



= 12.2. 



To fit the FRI observational data, it is necessary to intro- 
duce some luminosity evolution for these type of sources. 
However, Urry & Padovani (1995) report that the behav- 
ior of FRI radio galaxies is consistent with no evolution. 
Also, in the WJ scheme, the cosmological evolution of FRI 
does not explicitly appear, however, there are no funda- 
mental aspects in their model forbidding the evolution of 
FRI sources in general. The idea that FRI should evolve 
is confirmed by the recent work of Snellen & Best (2001), 
who on inspecting Hubble Deep Field images came across 
two sources of FRI morphology. This could mean that this 
kind of source is significantly more abundant at high red- 
shifts than at the present epoch. They suggest also that 
cosmological evolution of FRI sources could be similar to 
that of FRII (i.e. is radio luminosity dependent). We pro- 
pose that the evolution of FRI sources is similar to that 
of FRII, (Sect. 3.1, eq. 4). The evolution rate M, which 
depends on luminosity, has the form: 
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Table 3. Values of beaming and evolution parameters for 
models A and B. 



Parameter 


FRI 




FRII 






A 


B 


A 


B 


7 


10.3 (t) 


10.3 


20.0 (t) 


20.0 


RcM(*) 


0.010 (t) 


0.010 


0.004 (t) 


0.006 


Iog(Li(**)) 




26.50 


25.44 (t) 


24.50 


l0g(L2(**)) 






27.34 It) 


27.95 






5.0 


10.93 (i) 


12.20 


Zmax 


5.0 (t) 


5.0 


5.62 (t) 


5.62 



t the value is taken from Wall & Jackson (1997), and 
I the value is taken from Jackson & Wall (1999), 
{*) u = 5 GHz, (**) WHz-^sr-i 

M(L)=0 forL<Li, 
M(L)=Mmax forL>Li. 

We estimate that the small fraction of FRI with lumi- 
nosities Li5i > 10^®-'' (WHz~^sr~-^) evolve with cosmic 
time. This is supported by the fact that at L151 « 10^^'^ 
(WHz~^sr~^) in the low-frequency samples the fraction of 
objects with observed broad-line nuclei changes radically 
from ^ 0.4 at higher radio luminosities to ~ 0.1 at lower 
luminosities (Willott et al. 2000). 

The best-fit relativistic beaming parameters (Rc(i^) and 7) 
obtained are similar to those estimated by Wall & Jackson 
(1997). They are also in agreement with those derived 
from observations. The observed core- to-extended flux ra- 
tio of individual sources spans a wide range of values from 
~ 10~^ for the lobe-dominated sources (Morganti et al. 
1993) to ~ 10^ for the core- dominated sources (Murphy ct 
al. 1993). The observations show also that there is a large 
spread of Lorentz factor 7 values from 2 to 20 (Vermulen 
& Cohen 1994), while some numerical simulations give 
4 < 7 < 40 (Urry & Padovani 1995). However, the Lorentz 
factor value of FRII sources should be greater from those 
of FRI (Urry & Padovani 1995). 

Some recent papers (e.g. Caccianiga ct al. 2002) suggest 
the existence of another class of sources with compact 
morphologies that do not originate from FRI or FRII pop- 
ulations. Our investigation shows that there is no explicit 
need to introduce any additional population of compact 
powerful radio sources, apart from the beaming version 
of FRI and FRII parent populations. However, if the val- 
ues of beaming parameters given here are overestimated, 
then the model predictions for compact sources will be too 
small when compared to the data. 

5. Summary 

In this paper we have presented the observational log(n)- 
log(S) relation constructed for powerful extragalactic ra- 
dio sources separately for different morphological classes 
(i.e. FRII and FRI lobe-, and core-dominated). The 
counts are performed over three flux-density decades (from 
20 mJy to 20 Jy). We compare the observational source 
distributions with the predictions of the dual-population 



WJ model. From our analysis it can be inferred that: 

(i) the WJ dual-population evolutionary scheme enables 
predictions that are in good agreement with the observa- 
tional counts of powerful extragalactic radio sources of a 
given morphological type; 

(ii) we found a satisfactory set of physical parameter val- 
ues for relativistic beaming and cosmological evolution 
which in the latter case differ slightly from those estimated 
by WJ; 

(iii) to fit properly the observational data of FRI sources 
it is necessary to introduce an extra term of positive cos- 
mological evolution, i.e. in the proposed simple approach 
the most luminous FRI sources evolve in a similar way to 
the FRII sources. 
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Fig. 1. Normalized differential counts of sources of different morphological types at 1.4 GHz. Abscissa: log flux-density 
(Jy). Ordinate: log differential number of sources multiplied by S^'^[sr~^ Jy^-^]. Observational data are provided in the 
form of stepped curve with vertical error bars. The solid curve shows the counts predicted from model B and the dotted 
solid one - counts from model A. The respective panels give the counts for a) all sources in the sample, b) compact 
FRI and compact FRII sources, c) FRI lobe-dominated sources, and d) FRII lobe-dominated and CSS sources, (for 
details, see the text). 
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Fig. 2. Normalized differential counts of sources of different morphological types at 408 MHz. Abscissa: log flux- 
density (Jy). Ordinate: log differential number of sources multiplied by S^'^[sr~^Jy^-^]. Observational data are provided 
in tiie form of stepped curve with vertical error bars. The solid curve shows the counts predicted from model B and 
the dotted solid one - counts from model A. The respective panels give the counts for a) all sources in the sample, 
b) compact FRI and compact FRII sources, c) FRI lobe-dominated sources, and d) FRII lobe-dominated and CSS 
sources, (for details, see the text). 
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Fig. 3. Fraction of a) compact flat-spectrum sources, b) FRI lobe-dominated sources, and c) FRII lobe-dominated 
together with CSS sources in the sample at 1.4 GHz - left panels, and at 408 MHz - right panels. The solid curve 
shows the source ratio predicted from model B. 



